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Color codes are a class of topological quantum codes with a high error threshold and large set 
of transversal encoded gates, and are thus suitable for fault tolerant quantum computation in two- 
dimensional architectures. Recently, computationally efficient decoders for the color codes were 
proposed. We describe an alternate efficient iterative decoder for topological color codes, and apply 
it to the color code on hexagonal lattice embedded on a torus. In numerical simulations, we find an 
error threshold of 7.8% for independent dephasing and spin flip errors. 
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I. INTRODUCTION 

Topological quantum codes are well suited for fault- 
tolerant quantum computation [IH3] in two-dimensional 
qubit arrays constrained by short-range interaction. One 
of their topological features is that the number of encoded 
qubits depends only on the topology of the embedding 
surface but not on the code size. For such codes, the 
observables measured in order to identify errors are very 
close to local. They can be read out efficiently using local 
and nearest-neighbor quantum gates. 

The first family of topological quantum codes proposed 
were Kitaev's surface codes [4] . They combine many ad- 
vantages, namely a high error threshold [5 , very low 
weight of stabilizer operators, an efficient decoding al- 
gorithm [5j [6] , and the flexibility to implement encoded 
quantum gates in a topological fashion by code deforma- 
tions [7HTU] . A reason why one might want to go beyond 
surface codes is to reduce the operational cost of fault- 
tolerant quantum computation in two-dimensional qubit 
arrays. For surface codes, only a small set of gates can be 
realized transversally or by code deformation. All other 
gates invoke the procedure of magic state distillation [TT] 
which, although it does not dominate the overhead scal- 
ing, increases the operational cost. 

A further family of topological quantum codes are color 
codes [HI H3] • For them, a larger group of quantum gates 
can be implemented in a transversal or topological fash- 
ion, potentially reducing the operational overhead. For 
the color codes, the error-threshold turns out to be as 
high as for the surface code if the syndrome measure- 
ment is assumed perfect [14] . If error in the syndrome 
measurement is taken into account, surprisingly, color 
codes perform better than the surface codes [TS]. How- 
ever, for a full analysis in which the error-correction pro- 
cedure is broken down into gates, each with its own error, 
the error-threshold for the color codes seems to have the 
smaller value [16 . 

Error-correction for the color codes can be cast as 
a matching problem on a hypergraph [17] . for which — 
unlike the matching problem on graphs [3] that applies 
to the surface codes [5] — no polynomial time algorithm 
is known. Still, approximate polynomial time solutions 
can be found [TT] that lead to a good error threshold. A 



different approach is taken in [18] where the color code 
on a square-octagon lattice is mapped onto two copies of 
a surface code, and the surface codes are then decoded in 
a recursive fashion |19l somewhat analogous to concate- 
nated codes. Message passing is an integral part of the 
algorithm in [19] . 

In this paper, we present an alternate decoding algo- 
rithm for color codes, and apply it to the hexagonal lat- 
tice for which no threshold value has yet been reported. 
We find an error threshold of 7.8% for an error model 
with independent dephasing and spin flip errors. Like 
the algorithm in |19j . the present algorithm is recursive 
and employs message passing, but it decodes the color 
code directly without mapping to two copies of a surface 
code. 



II. BACKGROUND 

A. Topological color codes 

There are two main ingredients for constructing a color 
code: i) a trivalent graph ii) an embedding of the graph 
on a surface so that its faces are 3-colorable. The stabi- 
lizer generators of the code are defined as 

where / is a face of the embedding. 




FIG. 1. Topological color code on a torus 



In this paper we consider the color code on a hexag- 
onal lattice embedded on a torus. The same ideas can 
be applied to other color codes as well. We choose the 
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embedding shown in Fig. [T] The resulting code has the 
parameters [[n,4]], where n is the number of the vertices 
in the graph. The logical operators are associated with 
the nontrivial cycles on the torus, and carry a color label. 
The torus supports four encoded qubits. Our choice of 
hexagonal lattice and its embedding on the torus gives 
a code with the parameters [[18 • 4 m , 4, 2™ l+2 ]], where m 
can be any nonnegative integer. 



B. The dual lattice 

Although the algorithm can be designed for topological 
color codes on an arbitrary graph, the codes most likely 
to be used (especially in a fault tolerant setting) are those 
embedded in a regular, translation-invariant structure, 
i.e., a lattice . Three suitable lattices are known, namely 
the hexagonal lattice, square-octagon lattice (also called 
the truncated square tiling) and the truncated trihexag- 
onal lattice [5D]. We restrict our attention to codes on 
regular lattices, in particular to the hexagonal lattice. 
The color code on this lattice on a torus has to satisfy 
certain restrictions on its length in order to be embedded 
on the torus. Assuming that both the logical operators 
color code have the same weight, we have n = 2- (3x m ) 2 , 
where m is a nonnegative integer and x > 2. For the 
purposes of decoding it is helpful to represent the color 
code on the dual lattice. 




FIG. 2. Hexagonal color code on a dual lattice (opposite 
sides of the lattice are identified). In this picture each vertex 
corresponds to a check and each (triangular) face corresponds 
to a qubit. Further, logical errors appear as nontrivial cycles 
of faces. One such operator is highlighted by the shaded faces. 



into a smaller copy of itself. 

The first step in the decoding the color code is to di- 
vide the code into smaller code blocks each of which can 
be decoded using an efficient decoder. For instance, the 
[[72, 4, 8]] code shown in Fig. [2]can be divided as follows. 






FIG. 3. Dividing the lattice 



Each block can be decoded separately and then replaced 
by a fewer number of logical qubits to give rise to another 
color code of smaller size, which can then be decoded us- 
ing the same process. The checks in each block can be 
distinguished into three classes: i) those that are inte- 
rior to the block ii) those that are on the boundary and 
shared between two blocks and iii) those on the corner 
and shared between more than one block. 

In decoding the block individually, we use the checks 
interior to the block and those on the boundary but not 
those on the corners. These are used in the next level of 
decoding. The errors made in the lower levels of decoding 
are corrected in subsequent levels using these remaining 
checks. 

Recall that the qubit corresponds to a triangle. Thus, 
when we replace a set of qubits by another logical qubit, 
we expect it to emulate a qubit in the following sense: A 
single error on a qubit flips all the three checks on the 
corners of the qubit. Similarly, the rescaled qubit must 
be identified with a triangle whose corner checks can be 
flipped by a logical operator. The smallest such cell is as 
shown in Fig [4} 



III. THE DECODING ALGORITHM— HARD 
DECISION VERSION 

A. Rescaling 

The main technique at work in the present algorithm 
is a rescaling transformation in which the code is mapped 




Rescaled to 



FIG. 4. Rescaling a group of qubits. 



The shaded qubits in Fig. [4] correspond to a logical op- 
erator that emulates a bit flip on the rescaled qubit. This 
operator flips the checks on the corners of the rescaled 
qubit. The collection of these qubits then behaves as a 
single logical qubit with this logical operator. 




FIG. 5. Rescaled lattice, after applying the rescaling trans- 
formation of Fig. [4] to the lattice of Fig. [3] 



There are a restrictions on the scale factor for the 
rescaling. For instance, the cell 




is not legitimate for rescaling. 

The restriction arises for the following reason. In each 
cell, there are more error locations than syndrome bits. 
Therefore, each error can locally be corrected only up 
to an undetectable error X that remains for the next 
level. For the rescaling to work, the error must act as an 
encoded error on the rescaled face. That is, X either flips 
all corner syndromes or none. This observation leads to 

Lemma 1. Suppose that we have a hexagonal color code 
represented on its dual lattice. Then an equilateral tri- 
angular cell in the (dual) lattice of side n containing n 2 
qubits and n(n + l)/2 checks can be rescaled if and only 
if n ^ mod 3. 

Proof. Suppose that the cell has a length that is not a 
multiple of 3. The rescaled cell must have the following 
properties. There must exist an operator which acts like 
the bit flip operator on a single qubit, it must therefore be 
able to flip all the corner checks. The crucial observation 
that we need is that the two checks of the same color can 
be connected by error operators as shown in Fig. [6] 



Because of the 3-colorability of the color code, the 
checks on the three corners of the cell are of different 
color. Now assume that there is a single qubit error in the 
interior of the cell, (the checks could be on the boundary 
but the qubit itself is in the cell) . Then this error causes 
three checks (each of different color) in its corners to flip. 
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FIG. 6. 2-chains on the color code 




FIG. 7. Logical operator when n ^ mod 3 



Since all the cells are connected, it is possible to move 
these checks to the corners of the cell. The operator that 
achieves this acts as the logical operator for the cell. 

If on the other hand, the side is a multiple of 3, all 
the three corners of the triangle are of the same color. 
One 2-chain can be used to flip two corners. However, 
the third corner will either lead to a flip of one or more 
check in the interior of the cell if its support is restricted 
to the cell. So it is not possible for any error pattern to 
simultaneously flip all the three corners. □ 

The above result can be generalized to color codes that 
are not based on regular lattices. 



B. Splitting the syndrome 

The motivation behind the division of the lattice into 
smaller blocks is to decode each cell locally. Since the 
checks are shared between adjacent cells, this implicitly 
requires us to partition the measured syndrome between 
the adjacent cells so that we can proceed with the local 
decoding of cells. Each (non-corner) syndrome in the 
boundary of the cells must therefore be split between 
the adjacent cells. For concreteness let us consider the 
following situation shown in Fig. [8] Supposing a shared 
syndrome were measured to be s. Then the block on 
the top and bottom can contribute s* and s l such that 
s* + s l = s. 



For simplicity consider the following heuristic. Assume 
that we randomly choose the initial split of the check. 
Let p(e|sgSiS2) denote the probability of errors that are 
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FIG. 8. Splitting the shared syndromes 



consistent with the syndrome. Then the probability that 



1 is given by 



p(s t \sis 2 ) 



p(e\s t ,sis 2 ) 



p(e\s t = 1, sis 2 ) + p(e|s* = 0, sis 2 ) 



(2) 



Similarly, we can compute the probability of the lower 
cell syndrome s l as 



p{e\s l , s[s' 2 ) 



p(e\s< = ljsi^) +p(e\s l = 0,s[s' 2 ) 



(3) 



Next the cell on the top and bottom exchange messages 
so that their estimates for the split probabilities are con- 
sistent. We then enforce the condition that s — s* + s , 
which leads to an update of the p(s*) and p(s l ) as fol- 
lows. Under the assumption that s — 0, we compute 
p(s t \sis 2 s' 1 s 2 ) as 



p(s*|sis 2 )p(s'|sis 2 ) 



p^Ms'K'a) + (l-p(s*| ai a 2 ))(l -p(s l \s[s' 2 )) 
If s = 1, then we compute p(s t \sis 2 s' 1 s / 2 ) as 



p(s t |s lS2 )(l-p( S i | S ' 1 s 2 )) 



-KaVi-a)) + (1 -p(s*|s lS2 )Ms'K S2 ) 



(4) 



(5) 



After the conditional split probabilities have been ob- 
tained through Eqs. Q - ((Sj) we make a hard decision 
for the syndrome split variable on each cell boundary in 
the lattice. We either choose the most likely splitting, 
or make a biased random choice with the respective split 
probability as prior. Then, we repeat the entire proce- 
dure until the pattern of syndrome splits reaches a stable 
configuration. In practice, we iterate for a predetermined 
number of rounds. 



not accounted for. These checks will be used at the next 
level, to account for any errors in the estimates made by 
the decoding in the present stage. Before we can move to 
the next stage however, we need an error model. We next 
show how to derive this error model for the next stage, 
but first a few remarks on the complexity of decoding the 
basic cells. 

Because the cell is of a constant size as the length of 
the code increases there is no increase in the complexity 
of this decoder with the length. The overhead comes in 
the form of the number of cells which goes polynomially 
in the length of the code. For code of length n = 18- x 2n \ 
with a cell size of x, there are 0(n/x 2 ) cells. Therefore 
the decoding complexity is 0(n) at each level. The over- 
all complexity including the rescaling at each stage we 
obtain the complexity as 0(nlog 2 n). 



D. Rescaling the lattice 

For the iterative procedure to work, we need the er- 
ror model at each level of rescaling. This is obtained 
as follows. Based on the split probabilities we obtain 
a splitting of all the shared syndromes, which we shall 
henceforth call the reference splitting. This can be made 
randomly or simply choosing the most likely split for each 
shared syndrome. The latter is deterministic — the same 
set of split probabilities will return the same reference 
splitting. 

Assume that we have a certain splitting for the block 
being decoded. Then the block decoder returns an esti- 
mate for this syndrome and also a probability of error in 
the estimate. Now suppose that the syndrome is sqSiS 2 , 
then the rescaling algorithm makes the correction that is 
given by canonical estimate corresponding to So s i s 2 and 
returns the probability of error in the rescaled qubit as 



hard 
r rcscalcd,so s\si 



p(e + X\s s 1 s 2 ) 

p(e|s sis 2 ) + p(e + A|s S!S 2 ) ' 



(6) 



where e is the canonical estimate and e = e + X, the 
complementary error. We pass on this error probability 
as prior to the next level. 



C. Decoding the basic cell 

The basic cell is decoded using a complete lookup table 
which is an optimal decoder. Given the split probabilities 
after the previous step, we make a splitting for all the 
shared syndromes. Every cell in the code can now be 
decoded locally with the syndrome internal to the cell 
and the syndrome split on the boundary. The effect of 
this local error correction is that all syndromes in the 
interior and the boundaries of the cells are accounted 
for. Only the syndromes in the corners of the cells are 



IV. REFINEMENTS 
A. Soft splitting of syndrome 

The hard-decision heuristic for finding the configura- 
tion of syndrome splits described in Section [HI B| discards 
much available information. Namely, in each round of 
the iteration, split probabilities are being computed but 
immediately rounded to or 1. A better approach is 
the following soft decision version, in which a recursion 
relation is set up for the split probabilities. The split 
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probabilities are computed as 

p(e\s t ,s 1 s 2 )p(s 1 s 2 ) 



vis 1 ) = £ 

sis 2 

p(s 1 ) = E 



p(e\s l = l,sis 2 ) +p{e\s l = 0, sis 2 )' 

p(e|5 i ,s / 1 s / 2 )p(s / 1 s / 2 ) 

p{e\s l = l,s[s 2 )+p(e\s l =0,s' 1 s / 2 )' 



(7) 



(8) 



Once again, consistency with respect to the measured 
syndrome is enforced as in equations Q and ([5]). For 
instance, if the measured syndrome is s = we obtain 
the following update for the top split probability, 

pis t) _> P^)(l-P(s 1 )) (g) 

P{ ' p( S *)(l-p( S i )) + (l-p(s t ))p(s')' V) 

For the initial split of p{s l ) and p(s l ), we compute pi 
as probability of odd number of errors in e^eie^ and 
p r as the probability of an odd number of errors in the 
lower cell. Then we compute p(s*) and p(s l ) from pi and 
p r conditioned on the measured syndrome as in equa- 
tions Q and ([5]). 



B. Soft rescaling of the lattice 

As in the syndrome splitting case, the hard decision 
approach is not optimal. Since we do not know for sure 
which splitting is the correct one, it is preferable to use all 
the information available to perform the rescaling. This 
raises the problem of consistently combining the error 
probabilities from each splitting. 

The key idea in making this combination is the fact 
that a stabilizer generator has the effect of flipping the 
splitting on that check. This provides us with a means 
to relate canonically the errors across different splittings. 
Suppose that we have an estimate e for a splitting s. If 
we changed the splitting to some other splitting s, then 
the canonical estimate for this splitting is the same as 
the one obtained by applying half stabilizers that flip s 
to s. 

Let us elaborate on this in more detail. The choice for 
the estimate is follows. Let Si be the check operator on 
the zth check node (of the 2x2 triangular cell shown in 
Fig. [9j. Then the support of Si in the cell is given as 
Si| cc ii. For instance, the 0th check has support on qubits 
{0, 1,3}, they are shaded in the figure below. 




FIG. 9. Rescaling 



The crucial observation is that by applying So, we can 
flip the syndrome sq. So supposing we have a splitting 
so s i s 2 with estimate eoe^ea, to obtain the estimate for 
say S0S1S2, where sq = sq © 1 we can take the estimate 
of S0S1S2 and then flip the qubits on the support of So in 
the cell, thus we get e^e-x^^'i as the estimate for SqSiS 2 - 
So letting 0000 be the canonical estimate for the all zero 
syndrome we obtain the estimate for s§s\S2 = 100 as 
e§e\e 2 ez = 1101 and similarly for the rest. 

The rationale for this choice can be understood as fol- 
lows. We can view any splitting that differs from the 
reference splitting as having the same effect as the refer- 
ence splitting if it can be obtained from the estimate from 
the reference splitting by applying a stabilizer operator. 

Now we can give a soft decision based rule for updat- 
ing the rescaled error probabilities. Assume that a given 
cell has the reference splitting sqS\S2- If this splitting 
is correct, then we have the rescaled qubit's error prob- 
ability as given in ^. But supposing that the correct 
splitting was SQS1S2, then the error in the rescaled qubit 
will have to be conditioned on i) S0S1S2, and ii) the esti- 
mate for the cell and the correction were applied under 
the assumption that S0S1S2 was the correct splitting for 
the cell. 

Consider the stabilizer operator Sd 



Sd 



n a 



(10) 



i=0 

i:s;©s; = l 



Denote by = (Spleen) the restriction of Sd to the 
basic cell under consideration. The support of ed will 
determine if the erroneous splitting that was applied 
requires a correction on the rescaled qubit. Let e be 
the canonical error estimate associated with the split- 
ting S0S1S2 © SoSi§2- Then the error associated to S0S1S2 
which does not require a correction at the next level is 
given by e + e d . Then the probability of error in the 
rescaled qubit is given by 



hard 
/"rescaled , sq s i S2 



p{e d + e + X) 



p(e d + e + X) + p(e d + e) 



Prescaled ~ ^ ] -Prescaled,sosis 2 P( SoS l S 2) 



(11) 



(12) 



p(e d + e + X)p(s s 1 s 2 ) . . 
s^s 2 p(e d + e + X)+p(e d + e) 

Let us explicitly write down this rule for the case of 
2x2 cell. 



P 



p(1110) 
p(1110)+p(0000) 
p(0101) 



p(000) 



p(1001) 



p(1001)+p(0111) 
p(0010) 



TTT)P( °1)- 



^ p(0101)+p(1011)^ ulu ; + p(0010)+p(1100)^ Ui1 '' 
+ P(° 011 ) V (IQ()) + P(° 10 °) DflOlU 

+ pCooiiHpUioi)^ 1 ^ + p(oioo)+p(ioio)^^ lul -' + 

+ p( 1000 ) v (no) + p(1111) v(ni) 

+ p(1110)+p(0110)^ liU J + p(llll)+p(0001)^ lil J 

Note that the rescaled error probability is independent 
of the reference splitting. 
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C. Looking ahead at the corners 

For small cells, the algorithm has to slightly modified 
to account for certain malicious error patterns that would 
otherwise persist with appreciable probability at higher 
levels of the recursion. An example is 




Such patterns can be accounted for by modifying the 
qubit probabilities in the prior, by taking the corner syn- 
dromes into account "ahead of time". This will direct 
the decoding algorithm to correct these errors. 

Consider a corner syndrome. Each qubit in the check 
belongs to a different cell. 



\P3 













We update the probabilities of these qubits as follows. 
Denote by pi — p 6 , the probability of error in the qubit 
interior to the cell. Let s be the corner syndrome. Then 
the updated error probability of the qubit is given as 



Pe i-> 



PePi 



PePi + (1 — Pe)(l - Pi)' 



(15) 



where p e is the contribution of the qubits external to the 
cell 



\ -(-1)4 lit 1 - 2Pi)' 



(16) 



D. Belief propagation for prior update 

The performance of the algorithm can be improved by 
providing a more accurate estimate of the qubit error 
probabilities to the decoder. This can be achieved, as in 
the case of [19] by using a standard belief propagation 
decoder to update the prior estimates for the qubit error 
probabilities. Unlike the toric codes, the color codes have 
a smaller girth and a larger number of 4-cycles in their 
Tanner graph, therefore it is not recommended to use 
too many iterations in the belief propagation algorithm 
to update the priors. 

The complete algorithm for the decoder is now given 
in Algorithm [l] 



PS 
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5-HT 2 6-KT 2 7-KT 2 8-1CT 2 9 ■ 1(T 2 
Bit (Phase) Flip Error Rate 

FIG. 10. Simulation results. The logical error rate after error- 
correction vs the physical error probability of independent 
spin (phase) flip errors. Shown is the performance of [[18 • 
f, 4, 2 m+2 ]] codes for m = 2, 3, 4 and 5. 



Algorithm 1 Decoding topological color code 

Require: Syndrome, channel error rate. 
1: Update the initial error probability on the qubits. 
2: For each syndrome in the boundary compute initial split 
probabilities based on the qubits in the cell as per the 
following rule 



Send p or 1 — p to the neighboring cell depending on 
whether the shared syndrome is or 1. 
Enforce consistency on the split probabilities as per the 
following rule. 

PlPr 

VlVr + (l-pi)(l-p r ) 

where pi is the split probability as estimate by the cell 
and the p r is the message from the adjacent cell. 
Update the corner qubit probabilities based on the corner 
syndrome according to (16 1. Let pi be the probability 
that there is an error in the corner qubit and p e be the 
probability that there is an error in the qubits external 
to the cell. Then update the corner qubit probability by 
enforcing consistency on pi and p e as per ( |15[ l . 
Send messages to the neighbouring cells and estimate new 
split probabilities according to equations (JT]) , (JsJ and (J9j . 

Split the syndrome according to the split probabilities. 
Locally estimate the error and correct. 
Compute the error probabilities for the next level of re- 
cursion as in equation 1 13 1. 



Repeat until the level 0. 

Use a lookup table to estimate the final error. 
Propagate the error to lower levels to obtain the final 
estimate. 
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V. CONCLUSIONS 

We have presented a recursive decoding algorithm for 
the color code on the hexagonal lattice. It yields a thresh- 
old value of 7.8% for an error model of independent spin 
and phase flip errors, under the assumption that syn- 
drome extraction is perfect. The results of our numerical 



simulation are shown in Fig. 10 The recursion is based 
on a square cell of 8 qubits. 
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